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Abstract. We study the finite size effects on the thermai deconfinement phase transition 
from the hadronic gas phase to the QGP phase, using a simple thermodynamic model based 
on the coexistence of confined and deconfined phases in a finite volume. We first consider the 
bag model partition function for the QGP, then we include the necessary color singletness 
condition. By probing the behavior of some physical quantities on the hole range of temper- 
ature, it turns out that in a finite size system, all singularities occuring in these quantities 
in the thermodynamic limit (V — > oo) are smoothed out, and the phase transition is then 
rounded over a broadened critical region , with an additional shifting effect of the critical 
point induced by the exact color singletness requirement. An analysis of the finite size scaling 
behavior at criticality of some characteristic quantities allows us to determine the scaling 
critical exponents characterizing the thermal deconfinement phase transition. The obtained 
results are in good agreement with those predicted by other studies for a first-order phase 
transition. 



INTRODUCTION 



It is generally believed that at sufficiently high temperatures and/or densities a new phase of 
matter called the Quark Gluon Plasma (QGP) can be created. This is logically a consequence of the 
quark-parton level of the matter structure and of the dynamics of strong interactions described by 
the Quantum Chromo Dynamics (QCD) theory. Its existence, however, has been partly supported by 
Lattice QCD (LQCD) calculations and the Cosmological Standard Model predictions. Experimentally, 
the only way to study the QCD Deconfinement Phase Transition (DPT) occuring between the hadronic 
and QGP phases is to try to recreate, in Ultra Relativistic Heavy Ion Collisions (URHIC), conditions 
similar to those which occured in the early times of the universe (~ 10~ 6 ,s), immediately after the big 
bang. 

If the QGP phase is (will be) created in URHIC at SPS (LHC and RHIC) colliders, the volume 
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in which it takes (will take) place is certainly finite. Also, in LQCD studies, the scale of the lattice 
space volume is finite. This is our motivation to study the Finite Size Effects (FSE) on the expected 
DPT from a Hadronic Gas (HG) phase to a QGP phase. These effects are certainly important 
since statistical fluctuations in a finite volume can hinder a sharp transition between the two phases. 
Only in the Thermodynamic Limit (TL) are phase transitions characterized by some singularities: 
5— singularities for a discontinuous PT and power-law singularities for a continuous PT. In general, 
in finite systems, Finite Size Effects (FSE) lead to a rounding of these singularities and a possible 
mixing of phases. However, even in a such situation, it is possible to obtain informations on the 
critical behavior. Large but finite systems show a universal behavior called "Finite-Size Scaling" 
(FSS) allowing the put of all the physical systems undergoing a phase transition in a certain number 
of universality classes. The systems in a given universality class display the same critical behavior, 
meaning that certain dimensionless quantities have the same values for all these systems. Critical 
exponents are an example of these universal quantities. 

In the present work, we study the FSE on the thermal QCD DPT, using a simple thermodynamic 
model based on the coexistence of the HG and QGP phases in a finite volume, with low energy density 
phase and high energy density phase equations of state For the HG phase, we have calculated the 
Partition Function (PF) of a pionic gas as a first approximation, and for the QGP phase, we have 
calculated the PF of a free gas of quarks and gluons. This latter is first derived simply within the bag 
model, then we include the exact color-singletness condition, since the color confinement phenomenon 
requires that the QGP phase must have a neutral color charge. This condition is implemented into 
the quantum statistical description of the system using the Group Theoretical Projection Technique 
(GTPT) 2]. This Color Singlet Partition Function (CSPF) can not be exactly calculated and is 
usually calculated within the saddle point and gaussian approximations in the limit VT 3 >> 1 |3|]. It 
turned out that the use of the such obtained PF, within the Phase Coexistence Model (PCM), has 
as a consequence the absence of the DPT Q. This is due to the fact that the approximations used 
for the standard calculation of the CSPF, break down at VT 5 « 1, and this limit is attained using 
the PCM. We propose then a new method for calculating a suitable CSPF, containing an expansion 
in a power series of the QGP volume fraction, using the saddle point approximation but avoiding the 
problem arising at VT 3 « 1. The CSPF derived in this way allows us to accurately calculate physical 
quantities describing well the DPT at finite volumes, like the order parameter, the energy density and 
the entropy density, and probe the behavior of the PT at criticality. We show that in the TL, these 
thermodynamic quantities exhibit a discontinuity at a critical temperature T c (oo), which appears as 
a delta function singularity in the first derivatives, i. e. in the specific heat and the susceptibility. In 
a finite size system, the discontinuities disappear and the variations of these quantities are perfectly 
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smooth over the hole range of temperature. The delta singularities are smeared into finite peaks 
of width 5T(V), which is nothing but the broadening of the critical region, with the maxima of 
the peaks occuring at an effective transition temperature T C {V) shifted away from the true critical 
temperature T c (oo). This additional shifting effect of the critical temperature is induced by the exact 
color singletness requirement. An analysis of the FSS behavior at criticality of these maxima as well 
as of the width of the transition region and the shift of the effective transition temperature T C (V) 
relative to the true one T c (oo) allows us to determine the scaling critical exponents characterizing the 
thermal QCD DPT. 

2. THE QCD DECONFINEMENT PHASE TRANSITION IN A FINITE VOLUME 
2.1. Phase Coexistence Model and Partition Functions Calculation 

To study the QCD DPT in a finite volume, we consider a simple PCM used in .1], in which the 
mixed HG-QGP phase system has a finite volume: V = Vhg + Vqgp- The parameter f) representing 
the fraction of volume occupied by the HG: Vhg = §V, is then defined and can be considered as an 
order parameter for the DPT. The value: f) = 1 corresponds then to a total HG phase and when 
f) = 0, the system is completely in the QGP phase. Assuming non-interacting phases (separability 
of the energy spectra of the two phases), the total Partition Function (PF) of the system can be 
approximatively written as: 

Z{\ ] ) = Z QGP {\ ) )Z HG %. (1) 

The mean value of any thermodynamic quantity A(T, fj,, V) of the system, as defined in can then 
be calculated by: 

JA(\))Z(t))di) 

(A(T^,V))=°— , (2) 

;z(f,)df, 



where -A(fj) is the total thermodynamic quantity in the state f), given in the case of an extensive 
quantity by: 



Aft) = A HG {T, W) + A QGP (T, n, (1 — fy) V), 
and in the case of an intensive quantity, by: 

A(t)) = i)A HG (T, W) + (1 - tyA QG p(T, n, (1 - f)) V), 



(3) 



(4) 
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with Aqgp and Ahg the thermodynamic quantities relative to the individual QGP and HG phases, 
respectively. 

The partition functions of both HG and QGP phases are calculated, and their final expressions are 
given in the following. For a pionic gas, the PF is simply given by: 

Z HG = e£ T3v ™, (5) 

and for a QGP consisting of gluons and two flavors of massless quarks, with a chemical potential /i, 
within the bag model, the PF is given by 

ZeGP = exp((^ + ^ + ^^-|)^) , (6) 

where B is the bag constant, accounting for the real vacuum pressure exercised on the perturbative 
vacuum. 



2.2. Finite Size Effects on the DPT 



We'll study in the following the FSE on the thermal QCD DPT, at a vanishing chemical potential 
(/i = 0), and for this purpose, we'll examine the behavior of some thermodynamic quantities with 
temperature for varying volume. The quantities of interest are the order parameter, which is simply 
in this case the hadronic volume fraction fj (T,fj,,V), the entropy density s (T,n,V) and the energy 
density e (T,fi, V). Their mean values are calculated from eq.(2), and are given by the expressions: 

(S(T, M ,y))= fv{e fv_ 1) ( ? ) 

(8(T,IM,V)) = S QGP + ( 5HG ~ SQGP) (t}(T,lM,V)) 

(e (T, fj,, V)) = t QG P + {t HG - zqgp) <f> (T, V)}, (8) 



with: 



f — III _ 3771-M T 3 _ 2 T _ 1 , B 

J ~ I 30 90 1 P 1 1 2tt^ T ' T 



5qgp = Iff T 3 + 2^T; 5HG = ^T 3 (9) 

_ R i 37tt 2 rpj . 2t2 M 4 . „ _ IT 2 rni 



Fig.(l) shows the three-dimensional plots of the order parameter, the normalized energy density 
(e(T,V)) , , ,. , , . (s(T,V)) 

, and the normalized entropy density — ^ , vs temperature and system volume, at a 

vanishing chemical potential (/i = 0), using the common value B 1 ^ = 145MeV for the bag constant. 
The first-order character of the transition is showed by the step-like rise or sharp discontinuity of 
the illustrated quantities, approaching the TL, at a critical temperature T c (oo) ~ 104.5MeF, which 
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reflects the existence of a latent heat accompanying the PT. The quantities and ^j- are tradition- 
ally interpreted as a measure of the number of effective degrees of freedom; the temperature increase 
causes then a "melting" of the constituent degrees of freedom "frozen" in the hadronic state, making 
the energy and entropy densities attain their plasma values. In finite systems, the phase transition is 
rounded since the probability of presence of the QGP phase below T c , and that of the HG phase above 
T c are finite because of the considerable thermodynamical fluctuations. The transition region around 
the critical temperature T c (oo) is then broadened, acquiring a bigger width, smaller is the volume. 

3. THE QCD DPT IN A FINITE VOLUME INCLUDING THE COLOR-SINGLETNESS 

CONDITION 

3.1. The Color-Singlet Partition Function of the QGP 

The partition function for a color-singlet quark-gluon plasma contained in a volume Vqgp, at 
temperature T and quark chemical potential y, is determined using the Group Theoretical Projection 
Technique (GTPT) formulated by Turko and Redlich 2] and is given by Q{: 

Z QGP (T,f,,V QG p) = ^e- E ^ IIL J Jd^)d(fjMM)Z(T,tA,VQ GP] <p,i,), (10) 

— 7T —TT 

where M((p,ip) is the weight function (Haar measure) given by: 

MM) = (sin (I(V + §)) sin(f ) sin - f ))) 2 , (11) 
and Z the generating function defined by: 

Z(T, n, V QG p; <p, V>) = Tr [exp (-0 (h - fi(N q - Nq)) + i^T 3 + , (12) 



where (3 = — with the units chosen as: ks = h = c = 1, Hq is the free quark-gluon Hamiltonian, 



N q yNqj denotes the (anti-) quark number operator, and ^3 and Y$ are the color "isospin" and 
"hypercharge" operators respectively. 

The generating function Z (T, Vqgp, <p, ift) can be factorized into the quark contribution and the 
gluon contribution as: 

Z (T, n, Vqgp; <p, VO = Z quark (T, y,, Vqgp; <p, i>) Z g[uon (T, V QG p; <p, tp) , (13) 
where the quark contribution is given by: 



Zquark (T, fj,, Vq GP ; (p, tp) = exp 



q=r,g,b 



(14) 
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with q = r, b, g the color indices, dg = 2Nj counts the spin-isospin degeneracy of quarks and the 
angles a q being determined by the eigenvalues of the color charge operators in eq. (12): 



2 T 3' U S 



2 3' " b 



20 
' 3 ' 



and the gluon contribution is given by: 

Zgiuon (T, Vqgp; <p, tp) = exp 



7T 



7T 



%T 3 V QGP d G J2[~~5o + 

g=l V 

with d G = 2 the degeneracy factor of gluons and a g (g = 1, ...4) being: 

a.\ = a r — a g , «2 = oig — Qfe, 03 = at, — a r , = 
Let us write eq. (10) on the form: 



Z QGP (T^,V QGP ) = J_ e <vCi*»(*)-#) / / d(|) d (|)M( v3 ^)e^^)^^))^ 3 i 



(15) 



(16) 



(17) 



(18) 



with: 



fi _ 7T 2 21 

f' ~ 12^30' 



g(p^,^) = -(^-d Q + ^d G ) + -^- £ 



16 
15' 



TT^dQ 

12 2 



q=r,b,g 



-1 + 



("9 ~Kt)Y 



7T- 



19 9 2^ 



(a 9 - 7r) z 



5=1 



7T" 



(19) 



(20) 



go (y) being the maximum of g((p, tp, ^) for tp, tp € [— 7r, +7r], given as: 

„ I'M _ 7T 2 /21j 1 16 j \ 1 -f^h , 3fj? 3^_- 

fo It; - 12 \m d Q + T5"g?j + ^"l^r^ + 2i%f< 

and q the QGP volume fraction: q = 1 — fj. 

We evaluate the integral in eq. (18) by expanding M to leading order in ip and tp, and g to 
fourth order in ip and tp. Concerning the two terms appearing from g, we do an expansion of the part 
depending on the volume to second order for VT 3 ^> 1, while the term independent of the volume is 
expanded in a series on the form [J]: 

e »(^)q = ( g (^^)) J ' q .7 . (21) 

After some calculation, the final color-singlet partition function of the QGP, Zq G p, can then be 
expressed as p: 

4 e qV(T3 90 (#)-#; 



Z QG p(T,^,Fq) 



^2 /noc,n 
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12vr2 a 2 (^)yr3^ /piqP+1 

\ 1 / p=Q 



9- 2 (a(^)FT3) 4 [n=0 
where the ge neral form of the integral coefficients I nm is given by: 

n / 4 o/4 2/2\ m 

, / p 2tp ip tp \ 



dipdtpM^iip^) (ra!) _1 



1^ 



9 



+ 



22 



(23) 



with: a (^) 



if- ( + 1 J + f^G ) j M^ ' ) being the expansion of M to leading order in cp and 
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3.2. Finite Size Effects 

The mean values of the order parameter, the energy and entropy densities can also be computed in 
this case, using the obtained CSPF within the definition (2). Their expressions for two quark flavors 
are respectively given by: 

< h(T, V) > csc = 1 - V fLJ (24 ) 

< e{T, n, V) > csc = zqgp + e + (e ffG - egcp) < f>(T, n, V) > csc (25) 



< s(T, n, V) > csc = Sqgp + s + (5 H g ~ sqgp) < t)(T, V, V) > csc , (26) 

with: 

< '/m^)=((S+5o(^))T 3 -f) 

e = -12^ ;a = -^-^Iii(rr 3 a(^)) , 

ffo(r) ana ^ a (t) being taken for iV/- = 2. 

Fig. (2) illustrates as in the previous section the variations of the above quantities with temperature 
for different system sizes, but in this case including the necessary condition of color-singletness for 
the QGP. Additionally to the rounding of the transition in geometrically finite systems, the color 
singletness requirement causes a shift of the critical temperature to higher values, when the volume 
decreases. The internal degrees of freedom being gradually "frozen" with decreasing volume, the 
pressure of the color-singlet QGP, related to the energy density by the well known relation: Pqgp = 
3 ( e QGP — 4-B), is lower at a given temperature, and the mechanical Gibbs equilibrium between the 
two phases would then be reached for T C (V) > T c (oo) . For small size systems, the transition is then 
smoothed out over a range of temperature ST (V), around an effective transition temperature T c (V), 
shifted from the true transition temperature T c (oo). 

4. FINITE-SIZE SCALING ANALYSIS 
4.1. Finite-Size Scaling 

In the TL, phase transitions are characterized by the appearance of singularities in some second 
derivatives of the thermodynamic potential, such as the susceptibility x an d the specific heat c. 
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For a first order phase transition, the divergences are <5-function singularities, corresponding to the 
discontinuities in the first derivatives of the thermodynamic potential. For instance, in the case of the 
thermal QCD DPT, the ^-singularities appear in the susceptibility x {T, V) , defined to be the first 
derivative of the order parameter with respect to temperature, i. e., x {T, V) = gjT~ ~ ' an< ^ * n ^ ne 
specific heat density c(T,V) defined as: c(T,V) = — — — — . In finite volumes, these <5-functions 
are rounded into finite peaks over a range of temperature ST (V) as illustrated in Fig. (3), and the 
peaks occuring at an effective transition temperature T c (V) are shifted away from the true critical 
temperature T c (oo). The width of the transition region, the shift of the critical temperature and the 
maxima of the peaks of the specific heat and the susceptibility show a scaling behavior at criticality, 
characterized by scaling critical exponents as: 

ST (V) ~ V- 6t 

t t {V)=T c {V)-T c {^)~V- x t 
c™ ax (V) ~ V° T 

xT x {v) ~ v^t. 

The finite size scaling (FSS) analysis is used to recover the scaling critical exponents, and it has 
been shown 7] that for a first order phase transition, the finite size quantities ST (V), tt(V) scale as: 
V^ 1 , and the maxima c™ ax (V) and Xt &x (Y) scs ^ e as: V. The critical exponents 8t, At, 7t an d o-t 
are then all equal to unity. 



4.2. Numerical determination of the scaling critical exponents for the thermal QCD DPT 

4-2.1. The susceptibility scaling critical exponent J T 

The plots of the susceptibility are presented in Fig.(4-left) versus temperature for various sizes. The 
delta function singularity of the susceptibility occuring in the TL is smeared, in a finite volume, into 
a finite peak of width ST (V), with the maximum of the peak | Xt | max {V) occuring at the effective 
transition temperature T C (V). The plot of the maxima | Xt | max (Y) versus volume is illustrated in 
Fig.(4-right), and the linearity of the data with V can clearly be noted. A numerical parametrization 
with the power-law form: | Xt | max (Y) ~ V tT , gives the value of the susceptibility critical exponent: 
7 r = 0.99 ±0.04. Let us note that the associated error is a systematic one, resulting from the errors in 
the localization of the maxima. This systematic error will be estimated for all the critical exponents 
determined in the following. 
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4-2.2. The specific heat scaling critical exponent olt 

The variations with temperature of the specific heat density are presented for different sizes in 
Fig.(5-left), which shows the rounding of the delta function singularity of c(T,V) in finite systems 
into a finite peak of width 5T (V) and height c™ ax (V) . For decreasing volume, the width gets larger 
while the height of the peak decreases. The data of the maxima of the specific heat c™ ax (V) are 
fitted to the power-law form: c™ ax (V) ~ V° T in Fig.(5-right), and the obtained specific heat critical 
exponent is: ay = 0.99 ± 0.04. 



4-2.3. The shift scaling critical exponent \t 

For the study of the shift of the transition temperature tt{V) = T c (V) — T c (00), we need to have 
the finite size transition temperature T C (V). This latter is defined to be the temperature at which 
the rounded peaks of the susceptibility and the specific heat reach their maxima, and is found to 
be shifted away from the true transition temperature T c (oo). Fig. (6) illustrates the results for the 
shift of the transition temperature plotted versus inversed volume, and shows the linear character 
of the variations. The shift critical exponent obtained from a fit to the form: tt{V) ~ V~ Xt , is: 
A T = 1.0085 ±0.0009. 



4.2.4. The smearing scaling critical exponent 9t 

The width of the transition region can be defined by the gap: ST(V) = T2CV) — Ti(V) with T\(V) 

and T2{V) the temperatures at which the second derivative of the order parameter reaches its maxima, 

or in other terms the temperatures at which the third derivative of the order parameter vanishes, i. 

e.^{t>(T,V)) =0. 
dT ri(v),T 2 (v) 

Fig.(7-left) illustrates the variations of the second derivative of the order parameter with tem- 
perature for various sizes, and shows that the gap between the two extrema, which represents the 
broadening of the transition region, decreases with increasing volume. The results for the widths 
5T(V), plotted in Fig.(7-right) vs the inverse of the volume, were fitted to the power law form: 
8T (V) ~ V~° T , and the obtained smearing critical exponent is: 9t = 1.016 ± 0.007. 



5. CONCLUSION 



Our work has shown the influence of the finiteness of the system size on the behavior of thermo- 
dynamical quantities near criticality. The sharp transition observed in the thermodynamical limit, 
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signaled by discontinuities in the order parameter and in the energy and entropy densities at a criti- 
cal temperature T c (oo), is rounded off in finite volumes, and the variations of these thermodynamic 
quantities are perfectly smooth on the hole range of temperature. All discontinuities are rounded 
over a broadened critical region of width 5T (V) around the critical temperature, and the delta func- 
tion singularities corresponding to these discontinuities, appearing in the first derivatives of these 
quantities, i.e. in the susceptibility and specific heat density, are then rounded into finite peaks of 
widths 5T (V). The inclusion of the exact color-singletness requirement in the QGP partition function 
induces an additionnal shifting effect of the effective transition temperature T c (V) to higher values 
(T c (V) > T c (oo)). A FSS analysis of the behavior of the width of the transition region 5T (V), the 
shift of the effective transition temperature relative to the true one tt (V) = T c (V) — T c (oo), and 
the maxima of the rounded peaks of the susceptibility Xt^ x Y) an d the specific heat c™ ax (V) near 
criticality, shows their power-law variations with the volume characterized by the scaling critical ex- 
ponents 9t, At, olti It- Numerical results for these critical exponents have been obtained, and the 
associated systematic errors, resulting from the numerical method used for the localization of the 
maxima Xr^ (Y) an< ^ c T ax Y) an d °f Ti(V) an d T^Y), have been estimated. These numerical values 
are in good agreement with our analytical results: 0t = At = ar = 7t = 1 obtained in a parallel 



work 



These results are characteristic of the first order phase transition, as predicted by different 



FSS theoretical approaches 
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Fig. 1: Three dimensional plots of the mean values of (top) the order parameter, (middle) the 
energy density normalized by T 4 and (bottom) the entropy density normalized by T 3 , with 

temperature T and system size V, at jj, = 0. 
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Fig. 2: Plots of (bottom) order parameter, (middle) energy density normalized by T 4 and (top) 
entropy density normalized by T 3 versus temperature at /i = 0, for different system sizes including 

the color-singletness condition. 
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Fig. 3: Illustration of the critical behavior of the susceptibility x (T, V), the specific heat c (T, V) 
and the second derivative of the order parameter d\jdT . 




Fig. 4: (Left) Susceptibility x (T, V) as a function of temperature for varying volume, and (Right) 
linear fit of the results for the maxima of the susceptibility vs volume. 




5: Variations of (Left) the specific heat density with temperature for different system volumes 
and of (Right) the maxima of the specific heat density with the volume. 




Fig. 6: Plot of the shift of the critical temperature vs inversed volume. 




Fig. 7: (Left) Illustration of the second derivative of the order parameter (/i)" sc = d ^^Jf^ C3C vs 
temperature, taken at different volumes. (Right) Data of the width of the temperature region over 
which the transition is smeared, fitted to a power-law of the inverse volume. 



